Fractal structure of the effective action in (quasi-) planar models with long-range 

interactions 
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We derive the effective potential for composite fields in a class of (quasi-) planar models with long- 
range interactions. This class of models can be relevant for high temperature superconductors and 
graphite. The fractal structure of the effective potential is revealed and its physical interpretation is 
presented. It is argued that the multi-branched fractal structure of the potential reflects the presence 
of an infinite tower of excitonic bound states that occur as a result of the long-range interactions. 
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For a long time, there has been a great interest in low 
dimensional quantum field theories. Typical examples 
are the 1 + 1 dimensional Gross-Neveu model and the 
2+1 dimensional quantum electrodynamics (QED2+1). 
They have been used as very instructive toy models of 
relativistic systems with complicated dynamics. In addi- 
tion to that, it is quite remarkable that these same mod- 
els mimic some intrinsic properties of the dynamics of 
non- Fermi liquids. In particular, numerous recent stud- 
ies suggest that high-T c superconductors and graphite 
could be described well by a (quasi-) planar QED2+1 

mum urn iB. 

The difficulties in studying strongly correlated systems 
such as high-r c superconductors and graphite might be 
related in part to the presence of long range Coulomb like 
interactions between quasiparticles. This limits consid- 
erably the use of many powerful techniques that used to 
work well in models with short-range interactions. One 
of such tools is the celebrated Ginzburg-Landau (GL) ef- 
fective action. A method of the derivation of the GL ac- 
tion for composite operators in systems with short-range 
interactions has been known since the work of Gorkov 
0. Its generalization to systems with a spontaneous 
symemtry breaking driven by long-range interactions is 
far from being straightforward. In relativistic field the- 
ories, a method of the derivation of the effective action 
in models with long-range interactions (driven by gauge 
fields) was elaborated in Ref. 0. 

Recently, using the method of Ref. 0, the effective 
potential has been derived in the so-called reduced QED 
[3. In such a model, while quasiparticles are confined 
to a 2-dimensional plane, the electromagnetic interaction 
between them is three dimensional in nature. The model 
is relevant for a class of condensed matter systems that 
include highly oriented pyrolytic graphite among others 
In this letter, we will reveal a remarkable property 
of the effective potential in reduced QED, not considered 
in Ref. • It will be shown that the potential has a multi- 
branched fractal structure reflecting the presence of an 
infinite tower of excitonic bound states (a clear signature 



of long-range interactions). It is also important that, as 
was shown in Ref. :?], reduced QED is intimately con- 
nected with QED2+1, relevant for cuprates 0,1101 In 
fact, as will be discussed below, the qualitative features 
of the effective potential in reduced QED are common 
for a wide class of models with long-range interactions, 
in particular, for QED2+1. 

We start from a general definition of the effective ac- 
tion in a theory in which the spontaneous symmetry 
breaking phenomenon is driven by the local composite 
order parameter (0\if}i()\0) , where ip is the Dirac spinor 
of the quasiparticle field and {p is the (Dirac) conjugate 
spinor. Following the conventional way (see for example 
Ref. H3)i we introduce the generating functional W(J) 
for the Green functions of the corresponding composite 
field through the path integral: 



JW(J) 



DipDip exp 



J (x)ip{x)tp(x)\ 




(1) 



where J(x) is the source for composite field and L qp (x) 
is the lagrangian density of quasiparticles in the model 
at hand. Then, by definition, the effective action for 
the field a(x) = — (0|$^/>(a;)|0) is given by the Legendre 
transform of the generating functional W(J), 



r(cr) = W(J) - J d 3 xJ(x)a(x), 



(2) 



where the external source J(x) on the right hand side 
is expressed in terms of the field cr(x) by inverting the 
relation 



SW 
8J{x) 



<j{x). 



(3) 



The effective action L(<x) in Eq. provides a natural 
framework for describing the low energy dynamics in the 
model at hand. It is common to expand this action in 
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powers of space-time derivatives of the field a: 



I» = J £ 



-V(p) + ^Z'°'(<r)d lt <Td u <T + 



(4) 



where V(a) is the efective potential. The ellipsis denote 
higher derivative terms as well as contributions of the 
Nambu-Goldstone bosons if the latter are required by 
the Goldstone theorem. By making use the definition in 
Eqs. © and 0, we derive the following relation: 



ST 

5o{x) 



-J(x) 



(5) 



In the limit of a vanishing external source, this equation 
turns into an equation of motion for the composite field 
a(x). In a particular case of constant configurations, the 
equation reads dV/da = 0. 

By keeping the external source on the right hand side 
of Eq. JSJ nonzero but constant in space-time, we derive 
the following convenient representation for the effective 
potential: 



V{a) = -w(J) + Ja = da J {a), 



(0) 



where w{J) = W(J)/V 2 +i, and V 2 +i is the space-time 
volume. Thus, the crucial point for evaluating V is how 
the source J, playing here the role of the bare fermion 
gap, depends on the composite field a. It will be shown 
below, following the approach of Ref. , that the func- 
tion J(a) can be determined from a modified gap equa- 
tion, including the bare gap J. 

The model under consideration is the same as the 
model used in Ref. [a,l3 to study a metal-insulator phase 
transition in highly oriented pyrolytic graphite. We as- 
sume that the fermionic quasiparticles are confined to a 
flat 2-dimensional plane, while the electromagnetic field 
is free to propagate in the 3-dimensional bulk. The spa- 
tial coordinates on a plane are denoted by r = [x,y] and 
the orthogonal direction is labeled by the z coordinate. 

The lagrangian density of the electromagnetic field (in 
the bulk) is given by 

L em = -!- ( e Q E 2 - —B 2 ) - Aqp + -A ■ I (7) 
87T V Mo / c 

where £o is the dielectric constant, fi$ is the magnetic per- 
meability, Aq and A are the scalar and vector potentials, 
and electric and magnetic fields are E = — V Aq — \dtA, 
B = V x A. We note that the interaction terms, with 
the quasiparticle charge density p and current j, were ex- 
plicitly included in the lagrangian density in Eq. Q. In 
addition to that, the lagrangian density of quasiparticles 
themselves (living only on the plane) should be specified. 
It reads 



where vf is the Fermi velocity, ip(t,r) is a 4-component 
spinor, tJ) = V^7°, and the 4x4 Dirac 7-matrices furnish 
a reducible representation of the Clifford (Dirac) algebra 
in 2 + 1 dimensions. 

We consider the case when the fermion fields carry an 
additional, "flavor", index i — 1,2, ...,Nt. Then, the 
symmetry of the lagrangian JSJ is U{2Nf). Adding a bare 
quasiparticle gap (mass) term A^tptp into the lagrangian 
density (JSJ would reduce the U(2Nf) symmetry down to 
the U(Nf) x U(Nf). Therefore the dynamical generation 
of the gap would lead to the spontaneous breakdown of 
the U(2N f ) down to the U(N f ) x U(N f ). 

Proceeding as in Ref. |7j, we integrate out the bulk 
gauge bosons from the action. Then, up to relativistic 
corrections of order (vf/c) 2 , we are left with the following 
action of interacting planar quasiparticles: 



S qp ~ J dtd 2 rL a (t,f*)- i J dtd 2 r J dt'(fr'i()(t,r) 

x 7 V(*,rW - \r - r'|)<K* , ,f')7W,f'), (9) 

where the bare potential Uo(t, |r|) takes the following sim- 
ple form: 

rr / 1 —n \ e 2 6(t) f d 2 k .2tt e 2 8{t) 

C/b(t,M) = / T^exp zfc-r-W = 10 

|fc| eoln 



eo J {2nf 



Note, however, that the polarization effects may consid- 
erably modify this bare Coulomb potential. Thus, after 
taking this into account, the interaction is given by 



doj f d 2 k exp(— iut + ik ■ r) 



2n 



|fc| + n(u,,|fc|) 



(ii) 



where the polarization function H(uj, \k\) is proportional 
(with a factor of 2tt/sq) to the time component of the 
photon polarization tensor. 

The starting point in our derivation of the effective 
potential for the composite field a is the gap equation in 
a model with a nonzero constant external source J [see 
Eq. |JI||]. By noticing that the external source term Jijjijj 
enters the action in exactly the same way as the quasipar- 
ticle bare mass term, we easily derive the corresponding 
gap equation 0, 



J 



A = 



qdqA q K(p, q) 

V? 2 + (A q / VF y 

„2 



2(e v F + 7re 2 iV//4)' 



(12) 



where the approximate expression for the kernel K(p, q) 
is given by K(p, q) — 6{p — q)/p+0(q — p)/q- In the most 
important region of momenta |g| 3> A q /vp where the 
pairing dynamics dominates, the only role of the term 
(Aq/ui?) 2 in the denominator of the integrand on the 
Lq = VF^p(t,f) (i^°dt/vF — i^dx — T-l 2 dy) ip(t,r), (8) right-hand side of Eq. l|12|) is to provide a cutoff in the 
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infrared region. Therefore, one can drop this term, in- 
stead introducing the explicit infrared cutoff A/vf in the 
integral, where the gap A is defined from the condition 
A = A q=A / VF . This is the essence of the so called bi- 
furcation approximation. As a result, we arrive at the 
following equation: 



(13) 



where A ~ ir/a is the ultraviolet cutoff (a is a lattice size). 
The last integral equation is equivalent to the differential 
equation, 




P 2 A; + 2pA' p + XA p = 0, 



(14) 



with appropriate boundary conditions. A non-trivial so- 
lution to the gap equation, satisfying the infrared bound- 
ary condition, 



2 A ' I 



0. 



(15) 



P\p=A/v 

exists only for A > 1/4, and it takes the following form 

A 3/2 



A p = . 

sin dyjvpp 



sin 



V2 A 



(16) 



where v = V4A— T and 5 = arctani/. The critical value 
A c = 1/4 corresponds to a continuous quantum phase 
transition. 

In addition, the ultraviolet boundary condition, 



J= (A p + P A' p )\ p=/ 
produces the relation: 



J = 



A 3/2 



v vpA 
sin | — In — — 

sia(28)y/v^K V 2 A 



2S 



(17) 



(18) 



This relation determines J as a function of the gap A. If 
one supplies it with a similar representation for the com- 
posite field er(A), expression JfjJ) could be used directly 
to derive the effective potential, 



V(a(A)) 



(19) 



As follows from the definition, the expression for a as a 
function of A is given by the trace of the quasiparticle 
propagator (multiplied by the factor i). Thus, 



qdqA q 



Nf_ 

nv F Jo y/q 2 + (A q /v F y 



N 



f 



, P 2 A'(p) 
ttXvf 

N f A 3 / 2 VX 
7rAu^/ 2 sin(2<5) 



. , v v F A\ 
sm [_l n _l 



(20) 



Note that the second line of this equation follows from 
Eq. (|12|) after differentiating its both sides with respect 
to momentum p and substituting p = A. 

In order to get the most convenient representation for 
the potential, we introduce a new parameter Ao which 
denotes the dynamical gap in the case of a vanishing 
external source. It is determined from Eq. (|18JI at J = 0: 



Ao = Av p exp 



2ir~45 



(21) 



Therefore the critical coupling A c = 1/4 corresponds to 
a continuous phase transition of infinite order. 

By making use of the last equation, we obtain a 
rather convenient representation of the functions J(A) 
and cr(A): 



J(A) = 
a(A) = 



A 3/2 



sm(2S)y/v F A 
N f A 3 / 2 VX 
nXvp^ 2 sin(25) 



sm 



sin 



In -ir- 
2 A 



2 A 



(22) 



(23) 



Finally, by substituting these into Eq. (|19fl and perform- 
ing the integration on the right hand side, we arrive at 
the following parametric form of the effective potential: 



V(A) 



N f A 3 I 1 



2ttv 2 f 



2v 2 



1 — cos I v In 



A 



a(A) 



1 • ( i A 
— sm v In - — 

v \ A 

4N f A 3 / 2 VX 
tt(1 + v 2 )v% 2 L 

l-v 2 . (v, A 

sm — In — — 

2v \2 A 



(24) 



(v A 
cos — In — - 
\2 A 



(25) 



Around the global minimum, a = cr(A ), the potential 
is approximated as 



V(a) 



?r(l 



V 2 ) 2 V F 2 



48A/A 



lii - 



TO 



(26) 



This expression may suggest that there is nothing un- 
usual about the potential given by the parametric repre- 
sentation in Eqs. (|24|l and l|25l) . A closer look, however, 
reveals a rather rich, fractal, structure of V(a) with infi- 
nite number of branches near the origin a = (see Fig. ^ 
and discussion below). 

The extrema of the potential V(a) are determined by 
the equation dV/dcr — J(A) = 0. By solving this equa- 
tion, we obtain an infinite set of solutions for A: 



AL1 = A exp 

v 



71 = 0,1,. 



(27) 
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The corresponding values of a(A^ n ) are determined by 
Eq. J2HJ). It is easy to check that the second deriva- 
tive of the potential calculated at every extremal point 
in Eq. (J2ZJ is positive, d 2 V/da 2 > 0. Therefore, these 
extrema are local minima. The values of the potential at 
these minima cr(A|™^ n ) are calculated to be 

nA^) = "^(ALt) 3 <0. (28) 

We see that the solution with n — gives the largest value 
of the fermion gap and the lowest value of the potential. 
Therefore, it corresponds to the global minimum of V(o~), 
or in other words, to the most stable (ground) state in 
the model at hand. 




FIG. 1: The schematic fractal structure of the effective poten- 
tial V(a) that appears as the consequence of the long range 
interaction. 

It is natural to expect that the potential also has max- 
ima lying between those minima. However, the situation 
is more subtle: while, as a function of parameter A, the 
potential does have maxima, it does not have them as 
a function of a. Instead, the potential V(a) has turning 
points there, see Fig.^ Indeed, the first derivative of the 
potential (|24Jl with respect to the parameter A is zero at 
the following maxima: 

A \2 = A exp arctan y ' , (29) 

\ v v o — bis J 

(with n = 0, 1, . . .). At the same time, the first derivative 
of cr-field with respect to A is also zero at these points. 
As a result, the derivative of the potential with respect 
to cr, i.e., dV/da = j^/^, is nonzero there. 

From our analysis, we see that the potential V(a), de- 
fined parametrically in Eqs. Ij24(l and Ij25(l . is a multi- 
branched and multi-valued function of a. The location 
of the branching (turning) points are determined by the 
values of the gap parameter given in Eq. (|29|) . Each 



branch of the effective potential has a local minimum at 
a = a ^A|™^ . We also notice that, while the locations 

of the branching points converge to a = in the limit 
n — ► oo, the shape of all higher branches (after an overall 
scaling transformation) resembles the shape of the first 
branch. Therefore, the potential around the origin ex- 
hibits a fractal structure. This is shown schematically in 

Fig.m 

We will now demonstrate that the fractal structure of 
the potential reflects a rich spectrum of excitonic com- 
posites in this model. It is well known (see, for example, 
Ref. [10() that, because of the Ward identities for currents 
connected with spontaneously broken symmetries, the 
gap equation coincides with the Bethe-Salpeter equation 
for corresponding gapless Nambu-Goldstone composites 
(gapless excitons in the present case). Therefore, the in- 
finite number of the solutions a£™/ in Eq. I|27|) for the 
gap implies that there are gapless excitons in each of the 
vacua corresponding to different values of n. The gen- 
uine, stable, vacuum is of course that with n = 0. Let us 
show that the excitonic composites, which are gapless in 
the false vacua with n > 0, become gapped in the genuine 
vacuum with n = 0. Indeed, the transition from a false 
vacuum (with n > 0) to the genuine one corresponds to 
increasing the fermion gap, A^ n — > A^ n = A , with- 
out changing the dynamics. Therefore, as a result of this 
increase of the gap of their constituents, these excitonic 
composites should become gapped. Thus we conclude 
that the fractal structure of the potential reflects the 
presence of an infinite tower of gapped excitons that are 
radial excitations of the gapless Nambu-Goldstone com- 
posites. It is clear that this effect is intimately connected 
with the long-range interactions in this model. 

A confirmation of this point comes from studying the 
effective potential at finite temperature T and/or chem- 
ical potential /i. Because of the Debye screening, both 
T and /i lead to a dynamical infrared cutoff for inter- 
actions. Our studies show that the (dis-)appearance of 
the higher order branches is very sensitive to the values 
of temperature and the chemical potential. In particu- 
lar, all branches with n > no disappear at a temperature 
in the region < T < |A^, or at a chemi- 

cal potential in the region y^A^P^ < M < V^A^. 
This observation reflects the process of "melting" of those 
bound states whose binding energy is of order T (or jj,) 
or less. 

It is also clear that this picture should be quite gen- 
eral and valid for a wide class of models with long-range 
interactions. In particular, we have found that a simi- 
lar fractal structure takes place in the effective action of 
QED2+1 relevant for cuprates. 
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